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Most of the cosmological information extracted from the CMB has been ob- 
tained through the power spectrum, however there is much more to be learnt 
, from the statistical distribution of the temperature random field. We review 

some recent developments in the study of the Cosmic Microwave Background 
(CMB) anisotropies and present a description of the novel tools developed to 
^vq ' analyse the properties of the CMB anisotropies beyond the power spectrum. 1 
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1 Introduction 

r- -i 

C^) . The standard scenario of the universe includes an early inflationary phase dur- 

ing which the universe experienced a drastic expansion. As a consequence of 
this inflationary period the homogeneity, isotropy and flatness of the universe 
that we observe today can be understood as a natural outcome, without a fine 
tunning of the initial conditions. Besides, quantum fluctuations of the fields 
dominating the dynamics at that early phase constitute the seeds which, after 
Ci • gravitational growing, will form the large scale structure (LSS) of the universe 

at the present time (for a detailed description of inflation and its consequences 
see e.g. |141| ). The statistical distribution of the quantum fluctuations of the 
' field responsible for inflation, known as inflaton, in the vacuum state is Gaus- 

sian. In addition, the energy density fluctuations generated from them are also 
Gaussian since they are related by linear theory. Moreover, the anisotropies 
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This review was almost complete before the second release of WMAP data 
which was made publicly available on March 2006. The major improvement has been 
made in the polarization data which are now much better understood regarding both 
systematics and foregrounds. This has implied a significant change in the value of 
the optical depth (r ~ 0.09). Another difference with respect to the values of the 
cosmological parameters estimated from the first year data is the possible significant 
deviation of the spectral index from unity as derived from the second release of data. 
Appart from these changes the rest of the parameters as well as the deviations from 
Gaussianity ("anomalies") remain very much unaltered. 
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of the CMB are related to the energy density fluctuations by the linearized 
Einstein-Boltzmann equations, being therefore Gaussian distributed as well. 
This is a very important result since, for this particular statistical distribu- 
tion, all the properties of the temperature random field are included in the 
second order moment. In harmonic space this is called the radiation power 
spectrum, usually denoted by Ci. 

The cosmological parameters which characterize the universe leave their 
imprint in the radiation power spectrum. In the standard inflationary scenario 
we only need to measure this second order moment to obtain all the possible in- 
formation carried out by the CMB. This is in fact what many experiments have 
been doing with increasing precision in recent years. The BOOMERANG H2] 
and MAXIMA [HI] experiments onboard stratospheric balloons were among 
the first to establish the spatial flatness of the universe. Most notably, the 
NASA Wilkinson Microwave Anisotropy Probe (WMAP) has measured the 
Ci with high precision up to the multipole I « 600 |22I - This measurement 
has implied a high improvement in the determination of the cosmological pa- 
rameters, reducing the uncertainties below 10% |22()| . The uncertainties are 
further reduced when other cosmological data sets, like the LSS distribution 
of galaxies and SN la magnitude-redshift data, are combined with the CMB 
anisotropies \227\ . Moreover, the degeneracy among some parameters present 
when only a single data set is considered is broken. This concordance model 
has been further confirmed by the first detection of polarization fluctuations 
by DASI j!28j and the further determination of the polarization power spec- 
tra by WMAP [1221, DASI EEl, CBI Q]2] and BOOMERANG pSIITSlj . 
An independent piece of evidence comes from the cross-correlation between 
CMB anisotropies and fluctuations in the galaxy number density. Several au- 
thors have found positive cross-correlations between the WMAP data and the 
NVSS radio galaxy survey JS2], implying that the model of the universe is 
different from Einstein-de Sitter and, in the case of a flat universe, that the 
universe is dominated by a dark energy component [78J 11741 |SJ 12391 1161| . 

Measurements of the statistical distribution of the CMB temperature and 
polarization random field represent a unique probe of the early history of the 
universe. Thus, if the distribution is consistent with Gaussianity at a high 
precision then it would imply a nice confirmation of the standard inflation- 
ary model. On the contrary, if some significant deviations from normality are 
found then this result would contradict the standard single-field model and 
motivate the search for alternative scenarios of the early history of the uni- 
verse. There have already been many suggestions of the latter in the past 
based on many different physical phenomena. Most notably, non-standard in- 
flation models based on extensions of the standard one in many different ways 
-e.g. multi-field scenarios as the curvaton mechanism |145| . inhomogeneous 
reheating scenario, etc (see [31] for a review)- and cosmic defects including 
global and local strings, monopoles, domain walls and textures |240| . 

The structure of the paper is as follows. In Sec. El we summarize the most 
relevant physical effects which generate the CMB temperature and polariza- 
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tion anisotropies. Sec. |31 deals with the constraints imposed by recent CMB 
data sets on the cosmological parameters as well as other cosmological data 
sets combined with the CMB. In Sec. 0] we describe the properties of the 
CMB anisotropies beyond the power spectrum and discuss possible sources of 
non-Gaussianity, inhomogeneity and anisotropy in their distribution on the 
celestial sphere. In Sec. we discuss the methods developed for studies of the 
CMB temperature distribution and of its global isotropy whereas in Sec. El we 
review the analyses performed so far on CMB data sets measured by different 
experiments. Finally the conclusions are given in Sec. [7| 

2 Theory 

The CMB anisotropies are produced by different physical effects acting on 
this radiation before the last-scattering surface, when the universe was around 
380000 years old. These primary anisotropies contain very valuable informa- 
tion on the early history of the universe and can be observed today almost 
unaffected. After the temperature of the photons drop below some 3000 de- 
grees they are not able to ionize the hydrogen atoms and propagate freely. The 
physical effects producing the primary anisotropies at and after recombination 
can be summarised by the following equation [1501 12U4) : 

AT. , 4>Jn) „ f° d<j) , , . (AT. ,\ 

_( n )«^i_2 + 2 y -£dt + n-(v -v e ) + (—(n)) e (1) 

The gravitational redshift suffered by the photons in their travel from the last 
scattering surface to the observer is given by the first two terms in the r.h.s. of 
cq-ffl They are known as Sachs- Wolfe and late Integrated Sachs-Wolfe effect, 
respectively |202| . The velocity of the baryon-photon fluid at recombination 
generates a Doppler effect. The intrinsic temperature of the photons at recom- 
bination, represented by the fourth term in the r.h.s. of eq.^ a ls° contributes 
to the total anisotropy. Eq. ^ accounts for the anisotropies at recombination. 
Before recombination, changes in the gravitational potential, due to the im- 
perfect baryon-photon coupling, can also produce a gravitational redshift (the 
early Integrated Sachs-Wolfe effect). The accurate computation of all the con- 
tributions requires to solve the linearized coupled Einstein-Boltzmann equa- 
tions. Several codes have been developed to numerically solve those equations 
as a function of the cosmological parameters (about 10, see below), being some 
of them publicly available (e.g. CMBFAST [23], CAMB [THE])- A comparison 
among the results of different codes showed that the accuracy achieved is very 
good, reaching » 0.1% up to £ = 3000 |2T2| . 

As will be explained below, a very relevant statistical quantity to com- 
pute is the 2-point correlation function of the temperature anisotropies or 
equivalently the correlation of the spherical harmonic coefficients: 



(2) 



4 Enrique Martinez-Gonzalez 



where Cg is the anisotropy power spectrum and the a\ m are the coefficients of 
the spherical harmonic expansion 



The homogeneity and isotropy of the universe have been assumed in eq.0 i.e. 
the ai m are not correlated for different I or m. In Fig. ^ the scales at which 
the different effects dominate the power spectrum Cg are marked. As we can 
see, the gravitational effects dominate at the lower multipoles (large angular 
scales). At intermediate scales, 100 < I < 1000, the spectrum is dominated 
by several oscillations, usually called acoustic peaks. These peaks appear as a 
consequence of the balance between the gravitational force and the radiation 
pressure. At the smallest scales, I > 1000, the Cg are damped because of the 
width of recombination and the imperfections in the coupling of the photon- 
baryon fluid (Silk effect, |5T5] ). 



Fig. 1. CMB power spectrum Ce for the best fit model given in [22]. The gray 
band represents the cosmic variance. The spectrum has been computed with the 
CMBFAST code |2T5| . 

When the small scales in the matter distribution become nonlinear and 
their collapse give rise to the formation of the first stars and quasars they 
start to reionize the surrounding matter. This ionized matter can interact 
again with the microwave photons and produce new anisotropics. Secondary 
anisotropics can be either produced by the scattering with free electrons or 
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by the gravitational effect of the matter density evolution. The scattering can 
leave a very clear imprint when the photons happen to cross the core of rich 
galaxy clusters where the electron density and temperature are very high (of 
several 10~ 3 cm -3 and ~ 10 8 K, respectively). In this case the electrons inject 
energy to the photons through inverse Compton scattering, producing a dis- 
torsion of the black-body spectrum. This effect is called the Sunyaev-Zeldovich 
effect fm\ . 

The evolving gravitational wells in the large scale structure produce a gravita- 
tional redshift in the photons. This is known as the late Integrated Sachs-Wolfe 
effect |202j or the Rees-Sciama effect when the evolution is non- linear |200j . 
Besides, the trajectory of the photons is lensed by the same gravitational wells 
producing a noticeble effect on the CMB anisotropies at arcmin angular scales. 
One important issue in the observation of the CMB anisotropies is therefore 
to distinguish between primary and secondary anisotropies. In the case of the 
SZ effect, as the frequency dependence is different from the Planckian one, 
it is possible to separate the effect from the intrincsic CMB anisotropies by 
using multifrequency observations. However, for the lensing effect this is not 
possible and other properties, such as high-order correlations in the CMB 
temperature and polarization fields, have to be used (see e.g. |1()2I IM1 1214) 1. 




Fig. 2. Angular power spectra for temperature and polarization. The CVs are plot- 
ted for scalar perturbations on the left and for tensor perturbations on the right, 
where a tensor/scalar ratio of r = 0.01 has been chosen. The rest of the parameters 
have been fixed to the best fit model of [22] . All the spectra have been computed 
with the CMBFAST code. 

In the same manner that anisotropies in the temperature are expected in 
the gravitational instability scenario for structure formation, anisotropies in 
the polarization of the CMB are also expected. The physical effects giving rise 
to these anisotropies are, however, different. Linear polarization is generated 
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by Thompson scattering at the end of recombination when the growth of 
the mean free path of the photons allows temperature anisotropies to grow. 
In particular a quadrupole is formed in the reference frame of each electron 
producing polarization after the scattering (see, e.g., 03)- The expected level 
of polarization is however small w 5%. As shown by |254H113] two rotationally 
invariant quantities E, B can be constructed from the Stokes parameters Q, U. 
The different behaviour of E, B under parity transformations implies that only 
three spectra are needed to characterise CMB polarization: Cf \Cf , Cj E . 

In addition to the anisotropies generated by physical effects associated to 
the energy-matter density perturbations (scalar perturbations) , new anisotropies 
can be generated by a background of primordial gravitational waves (tensor 
perturbations) . This background is also a generic prediction of inflation. The 
primordial power spectra for both scalar and tensor perturbations are usually 
characterised by a scale-free law of the form: 

P s (k) = A s (k/k ) n %P t (k) = At(k/k ) n \k = OMMpc- 1 . (4) 

There is no observational evidence of the primordial background of gravi- 
tational waves yet. Upper limits have been imposed by the combination of 
WMAP and other high resolution CMB data with large scale structure data 
implying a tensor-to-scalar amplitude ratio r = A t /A s < 0.9 [22], see next 
sections. In any case, the amplitude of both temperature and polarization 
power spectra produced by tensor perturbations are believed to be several or- 
ders of magnitude below the ones corresponding to the scalar perturbations. 
However, the .B-mode can only be generated by gravitational waves and thus 
its detection represents a unique evidence of the existence of this primordial 
background. This is the reason why there is now so much interest in planning 
new more sensitive CMB experiments to detect polarization (see e.g. |30p. 
The detection of the .B-mode is expected to be much harder than the already 
detected .E-mode because of its intrinsically smaller amplitude and the rel- 
atively larger foreground emissions (expected from Galactic synchroton and 
thermal dust and extragalactic radio sources; for a recent discussion of these 
isues see 230 ). An additional complication comes from the lensing conver- 
sion of the .E-mode to the .B-mode 255 which is expected to dominate over 
the primordial B-mode at multipoles I > 100 |119j . Several methods have 
been developed to reconstruct the gravitational lensing potential from the E 
and B-mode polarization correlations and remove the lensing contamination 
|105l 11151 1211| . A more technical problem is the separation of the mixing of 
E and -B-modes for observations with partial sky coverage for which other 
methods have been proposed |139l I37S1 0B| • 

Outside the inflationary scenario there are other possibilities to generate 
B-modes. For instance, cosmic strings or primordial magnetic fields would 
generate a vector and a vector plus tensor components of metric perturbations, 
respectively, contributing to the .B-mode polarization (see e.g. |186l 1140] '). It 
is important to remark, however, that in the case of cosmic strings and other 
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topological defects their role as seeds for the large scale structure formation 
is already very much constrained by the observed Cg, see below. 

In Fig. |21 the different temperature and polarization power spectra corre- 
sponding to the primary anisotropies for the concordance model are plotted 
for the scalar and tensor perturbations assuming a value of r — 0.01. The 
power of the tensor perturbations is a few orders of magnitude below the cor- 
responding scalar power spectra for £ < 100 and decays strongly for larger 
multipoles. 




Many parameters are needed to characterize the standar cosmological 
model. Their variation changes the amplitude and shape of the temperature 
and polarization power spectra in many different ways (see Fig. [3] for their 
effect on the temperature power spectrum). The cosmological parameters can 
be classified depending on whether they characterize the background universe 
or the primordial power spectrum. 

The background Friedmann-Robertson- Walker universe and its matter and 
energy content are determined by the following parameters: physical bary- 
onic density, Wb — f2bh 2 ; physical matter density, w m — f2 m h 2 , where Q m is 
the matter density parameter including the contributions from baryons, cold 
dark matter and neutrinos, fl m = flb + ^CDM + Ov', physical neutrino density, 
w v = fi v h 2 ; dark energy equation of state parameter, w = pde/pde', dark 
energy density, Qde (in case w = — 1 the dark energy takes the form of a 
cosmological constant and its energy contribution is denoted by J?^) and the 
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Hubble constant, h = H /100km s _1 Mpc^ 1 . 

The scalar and tensor primordial power spectra are characterized by the fol- 
lowing parameters: amplitude of the primordial scalar power spectrum, A s , 
defined by P s (k) — A s (k/ko) ns , where fco = 0.05Mpc _1 ; scalar spectral in- 
dex, n s ; running index, a = dn s /d\nk, normally determined at the scale 
fco = 0.05Mpc _1 ; tcnsor-to-scalar ratio, r = A t /A s ; tensor spectral index, n t , 
which is normally assumed to be n t = —r/8 from the consistency relation of 
inflation (see e.g. |141| 1. 

There is an extra parameter that accounts for the reionization history of the 
universe: the optical depth, r = ot J t ° n e (t)dt, where ot is the Thompson 
cross-section and n e (t) is the electron number density as a function of time. 
Besides, there are two possible types of matter density fluctuations: adiabatic 
which preserve the entropy per particle and isocurvature which preserve the 
total energy density. The standard inflationary scenario predicts fluctuations 
of the adiabatic type. 

It is important to note that the manner in which some of the parameters 
enter in the calculation of the CV produce degeneracies. In particular, there is 
the well known geometrical degeneracy involving fl m , Qa and fl^, where flk 
is the curvature density parameter, = 1 — fi m — Qa (an example of this 
degeneracy can be found in Fig. 5 of |155| ). This fact makes clear the need to 
combine different cosmological data sets to break those degeneracies, as will 
be shown below. 

3 Cosmological constraints 

Several issues are important to bear in mind in the analysis of the CMB 
anisotropy data. First of all, the microwave sky is interpreted as a realization 
of the temperature anisotropy random field whose origin is in the inflationary 
era during the first moments in the history of the universe. The quantum 
fluctuations in the field dominating the dynamics at that epoch, called the 
inflaton field, generate energy perturbations which produce anisotropies in 
the CMB through the physical effects discussed in the previous section. The 
anisotropy field so generated is not ergodic on the sphere and so averages 
on the field do not coincide with averages on the sphere. This fact implies 
an intrinsic uncertainty in the statistical quantities extracted from the data, 
in particular in the power spectrum Ct, called the cosmic variance. For the 
standard model this uncertainty is easy to compute since the random field 
is expected to be Gaussian, ACt = C^/^/£ + 0.5 (this is valid for the whole 
sky, however when only a fraction of the sky f s k y is observed then the error is 

— 1/2 

increased by approximately f sky ) . The value of this uncertainty in plotted 
in Fig. ^ for the standard model. 

Second, at the millimeter wavelenghts where the intensity of the black-body 
spectrum with temperature ss 2.725K is maximum there are other astronom- 
ical foregrounds which are also bright. These are the Galactic foregrounds: 
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synchrotron, free-free and thermal dust, and also the extragalactic sources 
whose emission peaks in the radio or the infrared bands. The problem of sep- 
arating the foreground emission from the cosmic signal is crucial and requires 
multifrequency observations and a good understanding of the frequency and 
spatial properties of the foregrounds. In recent years a strong effort in this 
direction is being made by the CMB community, devising new instruments 
optimized to observe at frequencies around 100GHz where the CMB emission 
dominates over the foregrounds (for both temperature and polarization, see 
e.g. [230 ), and developing sophisticated methods of component separation 
(see e.g. [H7|E1 m this volume). In any case, there are still regions of the sky 
(near the Galactic plane and others dominated by extragalactic sources) where 
the CMB is unavoidably hampered by the foregrounds. These contaminated 
regions are usually masked and therefore removed from the data analyses. 

In the next subsections we will discuss the results obtained on the model 
of the universe based on the Ci. Following the standard model and some 
observational evidences, the anisotropies are assumed to be Gaussian and 
thus all the statistical information is included in the second-order moments. 
It is for this reason that most of the analyses performed to constrain the 
cosmological parameters are based on the Cg. Even the likelihood methods 
to derive those parameters are based on the Gaussian assumption (see for 
instance the review by 104 ). In Sec. and after we will also review the works 
made to test the Gaussian assumption. 

3.1 Cosmological parameters from the Cg 

In the last decade there has been a strong technological development in CMB 
devices allowing very sensitive measurements of the temperature of the mi- 
crowave sky. Due to the atmospheric emission at microwave frequencies the 
experiments have to be located at very dry places on Earth or onboard of bal- 
loons or satellites. In particular the space mission NASA WMAP has made 
possible to cover the whole sky from a priviledged position (the Lagrangian 
point 2 of the Sun-Earth system) , saved from contaminating emissions com- 
ing from the Earth, Moon and Sun. In tabled the values of the cosmological 
parameters determined with the first year data of that mission are shown. 
An interesting result that can be derived from this table is that the Einstein- 
de-Sitter model (i.e. flat geometry with null dark energy density) is many 
standard deviations away from the best fit model (asuming a flat universe 
and an optical depth r < 0.3) 220 . 

The resolution of the WMAP experiment allows a good precision in the 
determination of the Ct up to multipoles £ < 600. Other experiments based 
on the ground like ACBAR Q2E1, CBI [M], VSA and balloon-borne 
like BOOMERANG |112j are able to reach higher resolutions covering small 
patches on the sky with good sensitivity. The data obtained with these ex- 
periments therefore complement very well the WMAP ones. In Fig. 0] the 
temperature anisotropy results from WMAP and the other high resolution 
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Table 1. Cosmological parameters using only WMAP first year data. In the fit 
the universe is assumed to be spatially flat and the value of the optical depth is 
constrained to r < 0.3 (from |220p . The physical baryonic density parameter Wb 
is defined by wi, = Qbh 2 and similarly for the physical matter density w m which 
includes the contribution from all the matter species: baryons, cold dark matter 
(CDM) and neutrinos. 



Parameter Values (68% CL) 



w b 0.024±0.001 

w m 0.14±0.02 

h 0.72±0.05 
A s 0.9±0.1 

n s 0.99±0.04 



experiments as a function of the multipole are shown. We can see that the 
concordance model follows quite well the data up to f ~ 2000. For higher 
multipoles only ACBAR and CBI experiments have enough resolution and the 
data seem to indicate an excess of power as compared to the model prediction. 
This excess can be interpreted as a contribution from secondary anisotropics 
coming from the SZ effect |2E| or/and the emission of radio sources |2.31| . 
However, new multifrequency data at arcminutes resolution are necessary to 
confirm the excess and discriminate among possible causes. 
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Fig. 4. The temperature power spectrum Ce measured by WMAP, ACBAR, CBI, 
VSA and BOOMERANG, compared to the best fit model given by [22]. 



Since the recent detection of polarization by the DASI experiment 123 
several experiments have measured the TE cross-power spectrum Cj E (WMAP 
[T52] . CBI QH and BOOMERANG QH) and the EE one G f ( CBI HH2 
and BOOMERANG |168| b The sensitivity of the polarization data is not yet 
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Fig. 5. Polarization power spectra measured by DASI, WMAP, CBI and BOOB- 
MERANG (TE) and DASI, CBI and BOOMERANG (E). Also plotted is the best 
fit model given by |22|. 



comparable to the temperature one. However, it is already sufficient to con- 
firm the main features of the concordance model and, more specifically, the 
adiabatic nature of the primordial matter density fluctuations. Namely, that 
the peaks of the polarization power spectrum, Cf, are out of phase with the 
temperature Q ones. The TE cross-power spectrum, Cf E , and the E-mode 
polarization power spectrum, Cf , are shown in Fig. [S] 

The next challenge is the detection of the -B-mode polarization. As com- 
mented in the previous section this detection would unambiguosly indicate 
the existence of a background of primordial gravitational waves. Moreover, it 
is the only chance that we have to know about its possible existence in the 
next decade since experiments aimed to directly detect gravitational waves are 
still not sensitive enough. The amplitude of this background, A t , or equiva- 
lently the tensor-to-scalar ratio r, is proportional to the energy scale of in- 
flation. Whereas the all-sky ESA Planck space mission 2 |185| will be limittcd 
by instrument sensitivity, being able to detect values down to r a 0.05 230 , 
planned very sensitive ground-based experiments covering small patches of the 
sky and carrying arrays of 1000's of detectors like PolarBear |187j . Clover |T3rj| . 
BRAIN Q]2| or QUIET [T53] . are expected to be limitted by the ability to 
remove the foregrounds and are expected to reach values of r ss 0.01 in the 
best case. In order to significantly improve this limit we will have to wait for 



2 http: / /www.rssd.esa.int /Planck 
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the next generation of space missions with 10 3 - 10 4 detectors, now under dis- 
cussion by ESA, NASA and other space agencies [2461 I30| . The combination 
of complete sky coverage, very sensitive multifrequency observations and an 
optimal control of systematics make the space missions the ultimate experi- 
ments to go down to values r ~ 10~ 3 — 10~ 4 allowing to probe energy scales 
for inflation down to « 5 x 10 15 GeV [2301- 

3.2 Combination with other cosmological data sets 

The first observations showing the existence of a dark energy dominating the 
dynamics of the universe were those based on the luminosity distance-rcdshift 
diagram determined with supernovae SN la [1941 1181?] . Those results were, 
however, taken with certain caution because of the assumption made that the 
low and high redshift SN la had the same light curve behaviour (in addition, 
there were many other possible systematics which raised some concern). In 
any case, SN la data complements very well the CMB anisotropies since for 
the former the dependence of the model prediction on the dark matter and 
dark energy density parameters enters approximately as fl m — f^A, and as 
Q m + ft A for the latter (see Fig.EJ. 

The large galaxy surveys, 2-degree Field Galaxy Redshift Survey 2dFGRS 
[IHI] and Sloan Digital Sky Survey SDSS have recently reached w 200000 
measured redshifts on large fractions of the sky. A three-dimensional power 
spectrum P(k) of density perturbations is estimated with each data set. The 
cosmological parameters enter in the prediction of P(k) through the initial 
power spectrum (A s ,n s ) and the transfer function which linearly connects 
the initial and present spectra. Besides, an additional parameter, the bias 
b, is required to link the galaxy power spectrum to the matter one P(k). 
Although this parameter will in principle depend on the scale, however, it is 
found that at large scales b is scale-independent (see e.g. |227| ~). There are 
still other problematics that need to be corrected before using the data for 
accurate parameter determination such as the redshift-space distorsions due 
to galaxy peculiar velocities, survey geometry effects, ... Once those effects 
are corrected the galaxy redshift surveys play a key role in breaking the CMB 
degeneracies. For instance, the degeneracies of the model predictions on both 
data sets on the plane formed by the spectral shape parameter hf2 m and the 
baryon fraction Q^/ f2 m are almost orthogonal. 

Results on 11 free cosmological parameters combining CMB data (WMAP 
PI . CBI HSU, ACBAR H28) with the 2dFGRS galaxy redshift survey fTST] 
are given m table |U As can be seen from that table, the combination of differ- 
ent cosmological data sets has allowed for the first time an accuracy < 10% in 
the determination of most parameters. Results combining WMAP data with 
the SDSS galaxy redshift survey |22fi| have also produced similar values for 
the parameters |227| . 

In addition to CMB and galaxy surveys, other combinations including the 
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Table 2. Cosmological parameters from WMAP, CBI, ACBAR and 2dFGRS com- 
bined data (from |22|). 



Parameter Values (68% CL) 

w b 0.0224±0.0009 

,„ nl ot-+0.008 

w v < 0.0076 (95% CL) 

w < -0.78 (95% CL) 

Q DE 0.73±0.04 

h 0.7ltg : £ 

r 0.17±0.04 

A n oqq+0-086 

lis u - ooo -0.083 

n s 0.93±0.03 
-0.03lt°;™ 

r < 0.90 (95% CL) 



HST key project value for the Hubble parameter [HI], SN la magnitude- 
redshift data [19511229] . Lya forest power spectrum |f)()lll59) or abundancies of 
reach clusters of galaxies can help to improve the results |22UI 1211)1 11901 IllUj . 
Confidence contours in the plane (J? m , ^a) for the combination CMB+SN 
Ia+cluster abundancies are given in Fig. [B] The complementarity of the three 
data sets to break degeneracies is clearly shown. A similar effect in the plane 
(O m ,w) is obtained by combining CMB+2dFGRS+SN la (see Fig.Q. 

Recently, an independent piece of evidence that the universe is not 
Einstein-de Sitter (flat geometry with null dark energy density) has been 
found by cross-correlating the CMB map with galaxy survey maps. The same 
evolving gravitational potential wells which generate the large scale structure 
of the galaxy distribution also produce the gravitational redshift in the CMB 
photons at late times (the late ISW effect, see Sec.EJ. The amplitude and sign 
of the cross-correlation depends on three parameters floE-, and h. For a 
flat universe, as indicated by many different observations as discussed above, 
a positive signal will unambiguosly imply the presence of dark energy. This is 
the case found recently cross-correlating the WMAP map with different large 
scale structure surveys (the radiosource survey NVSS [S5], the X-ray survey 
HEAO-1 PHI, the optical SDSS [T], the near-infrared survey 2MASS gTg] ) 
|3T1 17511171 15112391. In mm three different methods were used to estimate the 
WMAP-NVSS cross-correlation: direct temperature anisotropy-galaxy num- 
ber density, cross-power spectrum and covariance wavelet coefficients. A clear 
positive signal was found in the three cases using a maximum likelihood analy- 
sis. The significance of a non-null dark energy reaches 3.5cr for the cross-power 
spectrum. This is also the maximum significance detection of the ISW effect 
up to date (see also |161j for a recent analysis using directional wavelets) . The 
reason for the different detection levels obtained with the three methods lies in 
the incomplete sky coverage available for the analysis. Otherwise, for a whole 
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Fig. 6. Confidence contours for the pair (J7 m ,J7yi) combining SN la, CMB and 
cluster density data (taken from the Supernova Cosmology Project). 



sky coverage the same detection level will be guaranteed by the maximum 
likelihood method since the likelihood is invariant under linear transforma- 
tions of the data. In Fig.|Hlwe plot the two maps together with the combined 
mask needed to avoid the regions contaminated by foregrounds in the WMAP 
map and unobserved regions of the NVSS survey. An analysis allowing varia- 
tions in the dark energy equation of state parameter w also shows a prefered 
value close to —1. The 2D confidence contours for the pair fioE,w and the 
marginalized likelihoods are plotted in Fig. [5] 

3.3 Summary of the main results 

From the previous discussion about the cosmological implications of the large 
amounts of data already collected on the CMB anisotropies, galaxy redshift 
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Fig. 7. Confidence contours for the pair (Sl m ,w) using either SN la, 2dFGRS and 
CMB or a combination of the three data sets (taken from the Supernova Cosmology 
Project). 



surveys, etc, the most important result is the convergence of all those cos- 
mological data sets towards the same model of the universe, the concordance 
model. Below we summarize the main characteristics of this model and other 
consequencies implied by the data: 

• The geometry of the universe is very close to flat. 

• The dynamics is dominated by a dark energy whose equation of state is 
almost that of a cosmological constant. 

• Most of the matter content is in the form of cold dark matter. 
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NVSS map 




Fig. 8. The radio sources NVSS and the combined WMAP maps including the 
common masks used in the cross-correlation analyses (taken from 239 ). 

• The large scale structure of the universe was seeded by quantum fluctua- 
tions in the very early universe that evolved via gravitational instability 
(the first evidence of this came from COBE-DMR |5T3] ). 

• The initial matter-energy density fluctuations were of the adiabatic type. 

• Topological defects did not play a dominant role in the structure formation 
of the universe. 

• Primordial gravitational waves do not appreciably contribute to neither the 
temperature nor the E-mode polarization anisotropics (up to the present 
sensitivities reached) . 

• The reionization happened at relatively early times, z > 10. 3 

• The initial fluctuations were close to a homogeneous and isotropic Gaus- 
sian random field (see however the next sections). 



Considering the second release of WMAP data, z is indeed around 10. 
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Fig. 9. The marginalized pdf for Ode (upper) and w (middle). Contours of the 
2D likelihood are given in the lower pannel (taken from 239 ). From left to right 
the three columns correspond to the results obtained with the WMAP-NVSS cross- 
power spectrum, covariance wavelet coefficients and temperature anisotropy-galaxy 
number density cross-correlation, respectively. 



4 Beyond the power spectrum: sources of 
non-Gaussianity 

As it was said before, single-field inflationary models predict Gaussian temper- 
ature fluctuations homogeneously and isotropicaly distributed on the sphere. 
In this case all the information is contained in the power spectrum or, equiv- 
alently, the two-point correlation function. In the previous sections all the 
discussion has been based on such statistical quantity. However, this pic- 
ture can be modified if the origin of the primordial fluctuations is different 
(e.g., the presence of topological deffects or non-standard inflation). If this 
is the case, the generated fluctuations will show non-Gaussian features (e.g., 
non-random phases) and the statistical information is not only contained in 
the second order moment. Even if the primordial fluctuations generated in 
the early universe are Gaussian, non-Gaussianities can appear after recom- 
bination when the non-linear matter density evolution takes place and new 
anisotropies are generated. These secondary anisotropies, as already discussed 
in Sec. sec:theory, can be produced by the gravitational lensing effect, the in- 
teraction of CMB photons with reionised matter or the gravitational effect 
of the non-linear evolution of the large scale structure. On the other hand, 
the universe is believed to be a Friedmann-Robertson- Walker (FEW) model 
with a simply-connected topology. Departures from this scenario, either its 
geometry or/and topology, will induce inhomogeneity or/and anisotropy in 
the statistical nature of the temperature distribution on the sphere (e.g. hot 
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spots, rings, multiple images). Systematics arising from instruments and data 
processing as well as residues left after foreground removal to get the intrin- 
sic CMB signal will inevitably introduce artefacts which will contaminate the 
primordial signal. It is therefore crucial to disentangle these artifacts from 
primary and secondary non-Gaussianity (see JS| f° r a recent study of the 
residues left by different component separation methods). 

4.1 Models of the early universe 

Even in the standard inflationary scenario small deviations from normality 
can appear as a consequence of second-order effects during the evolution of 
density fluctuations after inflation finishes (see |2J| for a detailed discussion) . 
The non-Gaussianity produced by those second-order effects must be therefore 
present in any inflationary model. In the case of non-standard inflation, as the 
curvaton scenario or inhomogeneous reheating, the primodial non-Gaussianity 
should be added to that term. The level of non-Gaussianity in the cosmolog- 
ical perturbations can be phenomenologically characterized by the non-linear 
parameter /jv l m the gravitational potential: 

$ = $L + f N L$L, (5) 

where is the gravitational potential at linear order. Although /jvl will in 
general depend on the scale, from a practical data analysis point of view it 
is useful to consider it as a constant [21]. Recently, a complete second-order 
calculation has allowed an estimate of an effective for some inflationary 
scenarios and considering only a limitted number of multipoles £ rn ax < 500 
due to computational cost |142| . The result is a value j e J[{lmax) — 4 with 
an increasing trend towards higher multipoles. This level of non-Gaussianity 
might be marginally detected with future experiments like Planck. 

Topological defects in the universe are expected to appear as a natural con- 
sequence of phase transitions in the early universe (for a detailed description 
see 12401 1. They are associated to symmetry breaking of fields responsable for 
the different particle interactions. According to their dimensionality they can 
be classified in the following types: monopoles, cosmic strings, domain walls 
and textures. We already know from the mismatch between the observed os- 
cillations of the Ci and the generic featureless spectra predicted by cosmic 
defects that their role in the generation of the cosmological fluctuations has 
to be sub-dominant. The contribution of both global and local defects to the 
C t is constrained to be < 10 - 20% [2J [7J E31 E52J EH] ■ However, their in- 
trinsically active role in seeding structure formation in the universe produces 
a characteristic non-Gaussian CMB temperature field. On large angular scales 
the superposition of many defects tends to Gaussianity by the central limit 
theorem. On scales smaller than the projected inter-defect separation non- 
Gaussianity can be seen as line-like discontinuities, in the case of strings, and 
as hot/cold spots in the case of monopoles and textures. Cosmic strings, the 
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best studied defects, have recently received renewed attention because they 
are produced in a wide class of string theory models where inflation arises 
from the collisions of brains. 

4.2 Geometry and topology 

The angular distribution of the CMB temperature and polarization fluctua- 
tions on the sky is very sensitive to the spacetime metric of the universe. On 
the largest scales it is assumed that the metric corresponds to the Friedmann- 
Robertson- Walker one. Deviations from this metric, even if small, can lead to 
noticeable signatures in the CMB anisotropy. In particular for homogeneous 
and anisotropic models (classified as Bianchi models), it is well known that 
global shear or rotation can produce a spiral pattern or hot spots on the 
CMB sky |19l ll(J8j . Since those signatures are on large angular scales, strong 
limits on those models were already derived from the COBE-DMR experi- 
ment [1211 11521 I35| . More recently, asymmetries/non-Gaussianities found in 
the WMAP first-year data have been nicely fitted with a Bianchi VII/i model 
108 . However, even considering extensions of the Bianchi VII/! models in 
which a Qa term is included, the values of the parameters which fit the large 
scale asymmetries are ruled out by current observations at high significance 

Due to the local character of the General Theory of Relativity the global 
topology of the universe is not theoretically constrained (for different topolo- 
gies discussed in the literature see e.g. 129, 244, 23*51 1137| ). Non-trivial topolo- 
gies leave imprints in the CMB anisotropy at large angular scales which can be 
observed as anisotropic patterns, matched circles or, more generically, devia- 
tions from a Gaussian random field (see e.g. [58111061 1197| ). The observability 
of those signatures depends on the topology scale as compared to the horizon 
scale [206] . Detailed CMB simulations in non-trivial topologies are required 
to extract the maximum information from the data (see |193l 110(7) ). Preci- 
sion all-sky surveys are necessary to obtain strong limits on global topology. 
The COBE-DMR data already constrained several topological models (see 
e.g. [6"51 IT551 More recently, the WMAP data were also analysed to fur- 

ther constrain topology finding no evidence of non-trivial topologies |127M59| . 
Moreover, theoretical limits are derived for the size of the topologies that can 
be detected |127| . Besides, [1441 11991 II lj have claimed a dodecahedral topol- 
ogy for the universe to explain the anomalies found at the lowest multipolcs of 
the WMAP data. The next all-sky experiment, the Planck mission, will obtain 
more sensitive data in a wider frequency range allowing a better separation of 
the CMB from Galactic emissions and therefore a better determination of the 
low-order multipoles. It will thus help to clarify the present situation and to 
determine whether the present anomalies are due to systematics, foreground 
residues or are indeed true anisotropies of the CMB. 
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4.3 Secondary Anisotropies 

As already commented in Sec. |21 the temperature anisotropies of the CMB 
are usually divided in primary and secondary. The former are generated prior 
to recombination and are directly related to the initial density fluctuations. 
The latter are generated after matter-radiation decoupling and arise from 
the interaction of the CMB photons with the matter. These interactions can 
be of gravitational type (e.g., Rees-Sciama effect, |2UUII15U| V or of scattering 
type when the matter is ionised (e.g., Sunyaev-Zel'dovich (SZ) effect of galaxy 
clusters |221) or Ostriker-Vishniac effect for a homogeneous re-ionisation [1761 
1241) or, on the contrary, inhomogeneous re-ionisation of the universe [3J 1881 
1118) ~). Secondary fluctuations will in general induce non-Gaussian signatures 
that will add to the intrinsic statistical properties of the CMB anisotropies. 
Their contributions to the bispectrum (third order moment in harmonic space, 
see next section) during reionization have been studied in [54) . More recently, 
it has been shown that the imprint of the Ostriker-Vishniac effect can be 
characterized by an undetectable bispectrum and a significant trispectrum 
(fourth order moment in harmonic space) 0O| . The Rees-Sciama effect due to 
non-linear evolution of matter is expected to produce three-point correlations 
much below the cosmic variance [164) . 

The SZ effect due to galaxy clusters is one of the most important sources 
of secondary anisotropies. Two types of SZ effect can be distinguished: the 
thermal effect, induced by the CMB photon scattering off free electrons in 
the hot intra-cluster gas, and the kinetic SZ effect due to Doppler shift of the 
photons when the clusters move with respect to the CMB rest frame. Owing 
to its peculiar spectral signature the thermal effect will be separated to a 
good accuracy and thus removed from the cosmological signal. However, the 
contribution from kinetic SZ effect, which is spectrally indistinguishable from 
the primary anisotropies, will not be easily subtracted, and will thus remain 
a potentially significant non-Gaussian foreground contribution. The SZ effect 
of galaxy clusters is by nature a non-Gaussian process. Several studies have 
aimed at characterising the non-Gaussian signature, either in wavelet space 
(see e.g. 0]) or in the real one ([SH1E53 ), in order to address its detectability 
and extract the maximum information from it. 

On its way from the last-scattering surface to us, the CMB radiation passes 
mass inhomogeneities which deflect its paths by the gravitational lensing ef- 
fect. This process can approximately be described as a random walk of light 
through a continuous field of mass inhomogeneities, and thus leads to a diffu- 
sion process. The diffusion tends to broaden structures in the CMB on angular 
scales smaller than ~ 10' (see for a review). 

As long as the density inhomogeneities can be described as a Gaussian 
random field, the light is slightly deflected and its distribution is a Gaussian 
random process. A purely Gaussian CMB would therefore remain approxi- 
mately Gaussian despite the deflection at least for I < 1000. However, at 
smaller scales non-Gaussianity is imposed on an originally Gaussian CMB 
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by producing correlations between temperature fluctuations at large scales 
(through the late ISW effect) and light deflections which generate noticeble 
effects at several arcmin resolution |1U2I 1841 I214| . Those correlations induce 
high-order correlations in the CMB temperature and polarization fields. In 
particular, a bispectrum is generated whose signal to cosmic variance noise is 
higher for polarization and polarization-temperature bispectra than for those 
involving the temperature alone |1()2| . 

Gravitational lensing by galaxy clusters can also lead to non-Gaussian fea- 
tures in the CMB through their strong lensing effect. While their origin is not 
physically distinct from lensing by large-scale structures, galaxy clusters can 
be individually identified through their thermal Sunyaev-Zel'dovich effect and, 
for the brightest cases, extracted from CMB maps. Methods to remove non- 
Gaussianity from CMB maps therefore have to distinguish between lensing by 
clusters and by large-scale structures. 

5 Methods to detect and characterize deviations from 
Gaussianity, standard geometry and trivial topology 

There are infinite ways in which a random field can departure from Gaussian- 
ity. For instance, the cumulants above order 2 should be zero for a Gaussian 
distribution and thus a non-null value taken by any of them would represent a 
deviation from normality. In Fig. HOl it is shown a Gaussian map correspond- 
ing to the concordance model together with other non-Gaussian ones with the 
same power spectrum but with small amounts of skewness or kurtosis. Hence, 
there is not a unique way to detect and characterize deviations from Gaus- 
sianity. Depending on the kind of features that are investigated some specific 
methods will prove to be more efficient than others. A variety of methods have 
been already proposed to search for non-Gaussianity in the temperature or 
polarization maps. Acting in different spaces (real, harmonic, wavelet, eigen- 
mode, ...) they are able to extract relevant information which is otherwise 
hidden in the fluctuation maps. Improved methods for the combined temper- 
ature and polarization maps are expected to be developed in the coming years 
motivated by the large amounts of new data expected to be collected by the 
new experiments (specially the all-sky WMAP and Planck satellites). 

An important difference among the methods is wether they are guided or 
blind. The former are intended to test specific models of non-Gaussianity. The 
latter do not assume any form of non-Gaussianity. Guided methods would be 
more powerful in constraining non-Gaussianity than the blind ones at the cost 
of being model-dependent. Both methods have been used in the literature with 
different results: while guided methods have imposed upper limits on the non- 
linear parameter f^L some blind methods have obtained significant deviations 
of Gaussianity (see Sec. 16.21) . 

Below we summarize different methods which can be used to detect and 
characterize non-Gaussianity in temperature/polarization maps. 
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Fig. 10. Map simulations for Planck LFI 33 GHz (7V sid c = 256, FWHM = 33 ar- 
cmin) . Non-Gaussian features have been introduced in the map using the Edgeworth 
expansion (from |153p . 
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Methods based on real space are very useful for picking out statistical 
features. However, they may be insensitive to any localized features. Standard 
statistics like moments and cumulants, obtained from the moment generating 
function and its logarithm or cumulant generating function, are included in the 
1-point probability distribution function (pdf). Any non-zero amplitude for 
any cumulant of order higher than 2 would be an indication of the presence of 
non-Gaussianity. It is more common to use the skewness and kurtosis, the third 
and fourth standardized moments, which are defined by the ratio of the 3 and 
4-order cumulants nomalized to the square root of the 2-order cumulant to the 
power of the corresponding order, respectively. All the statistical information 
is, however, included in the n-point pdf. For the Gaussian case, it is sufficient to 
specify the 2-point correlations. Any positive detection of reduced correlations 
above order 2 would imply deviations from Gaussianity in the data. Tests of 
normality based on the Edgeworth expansion for estimating the skewness of 
the COBE maps have been developed by [53]. The Edgeworth expansion allows 
to describe perturbations of the Gaussian distribution in terms of a series of 
cumulants (see |114j ). An alternative expansion derived from the Hilbert space 
of a linear harmonic oscillator, producing a proper normalized distribution 
under any truncation, has been proposed by |19ri| . Other approaches based 
on goodness-of-fit tests of normality have recently been developed and applied 
to the MAXIMA data |4*2"] . 

A different approach to test Gaussianity is based on scalar quantities on 
the sphere. These quantities are constructed with the first and second covari- 
ant derivatives of the temperature field. A well known example is the local 
curvature (see [73 EH) but many others can be defined as the modulus of 
the gradient, shape index, determinant of the Hessian matrix, ... (see |166j 
for a detailed study of those quantities). These scalar quantities have been 
recently shown to be good detectors of small deviations from Gaussianity as 
given by the Edgeworth expansion |lri7) . Analyses based on statistics of ex- 
trema and excursion sets (e.g., number, correlations, eccentricities) provide 
very robust results since extrema are often many sigmas above the noise level. 
Those quantities can be computed (semi-)analitycally for a Gaussian field 
as was shown in the pioneering work [2?J- Also, it has been shown that the 
Gaussian curvature of peaks is an efficient discriminator of non-Gaussianity 
|lfi| . The extrema correlation function for a Gaussian field is accurately com- 
puted in O E7J whereas ^3] computed the correlation of excursion sets for 
Gaussian and some non-Gaussian temperature distributions. Minkowski func- 
tional describe the morphology of the excursion sets which are defined by 
thresholding the temperature field. They are three for the 2D sphere: total 
contour of the excursion set, total area and the total curvature (or the genus). 
Analytical expresions can be obtained for the Gaussian field. These morpho- 
logical descriptors have been implemented on the HEALPix pixelisation 85 
of the sphere by |242j . 



24 



Enrique Martinez-Gonzalez 



Alternative methods based on multifractal analysis and surface roughness, 
widely used in many branches of physics, have also been proposed to detect 
non-Gaussianity in CMB maps |7UII165| . 

Spherical harmonics space statistics are attractive as the predictions for 
Gaussian fields can in some cases be computed analytically, and their covari- 
ance properties may also be computed in principle. In practice, we may not 
be able to benefit from these advantages due to the complexities and size of 
the experimental data sets. As for the other statistics, sensitivity to artefacts 
such as those arising from the observing strategy and, in particular, incom- 
plete foreground subtraction, need to be assessed. The bispectrum, £?^ 2 ^ 3 , 
is the rotationally invariant third-order moment of spherical harmonic coeffi- 
cients and has been very often used to test different data sets. It is given by 
the following expression (see HT\ for a derivation of this and higher order 
moments) : 



where (...) is the Wigner-3j symbol. The bispectrum must satisfy the selection 
rules, requiring that £\ + £2 + £3 be even, the triangular rule \£i — £2] < £3 < 
£\ + £2 and that 7713 = m\ + 777,2- For simplicity reasons often the case that 
all multipoles are equal, £\ — I2 — £3, has been considered (see e.g. [ZZj for 
an application to COBE-DMR). Inter-scale correlations of the form Bi-ue+i 
have also been used (e.g. |147j ). They are generically expected in models 
with scaling seeds like cosmic defects. Other components with more distant 
multipoles may be considered but they are often dominated by noise. A fast 
algorithm to compute the bispectrum from HEALPix-formatted data has been 
developed by and applied to the HEALPix-formatted COBE-DMR data. 
|126| have recently implemented an algorithm to compute the trispectrum 
(fourth order moment in harmonic space) in the HEALPix pixelization and 
have applied it to the COBE-DMR data. 

Novel methods have been recently developed based on other spaces where 
the data are transformed by either filters, wavelets or eigenvectors. Since the 
operations involved are linear the new coefficients should form a set of Gaus- 
sian variables if the underlying map is Gaussian. Wavelets are compensated 
filters which allow one to extract information which is localized in both real 
and harmonic space and for this reason they can be more sensitive than clas- 
sical methods (for a detailed review of the properties of wavelets see 
In recent years they have been applied to the non-Gaussian analysis of the 
CMB. Many of the statistics and techniques used on real or harmonic space 
can be translated and applied to wavelet space (e.g., moments, scale-scale 
correlations, 1-point pdf). On small patches of the sky several planar wavelet 
families have been successfully applied to detect simulated cosmic strings em- 
bedded in a Gaussian signal ^H] and search for discontinuities due to sec- 
ondary anisotropies |3|. The Haar wavelet (the simplest wavelet that one can 
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construct) |228l and the Mexican Hat wavelet (given by the Laplacian of 
a Gaussian) [1531 |4"T] have been now implemented on the sphere (see Fig. ^2 
for an application of the latter wavelet to COBE). |153 | has shown that the 
spherical Mexican Hat wavelet can be more sensitive than the spherical Haar 
wavelet for detecting types of non-Gaussianity expected in some non-standard 
scenarios. Moreover, the former wavelet has been found to be more sensitive 
than the bispectrum for constraining the nonlinear coupling parameter ap- 
pearing in non-standard inflation using the COBE/DMR data |48| . Recently, 
directional wavelets having a prefered direction have been implemented on 
the sphere -e.g. elliptical Mexican Hat and Morlet- and applied to the CMB 
analyses |160| . The computational cost of exploring all possible directions can 
be greatly reduced by using steerable wavelets where the coefficients at any 
direction can be obtained as linear combinations of the ones of a basis with a 
limitted number of elements [2471 1248j . 

Filters can also help to extract non-Gaussian features from a given data 
set. This can be done in different ways, e.g. by removing the non-cosmological 
information (Wiener), equilizing the CMB signal (signal-whitening) or the 
noise (noise-whitening) 249 or taking field derivatives to study discontinuities 
(for an application to cosmic strings see 86 ) . In the noise- whitening case the 
corresponding signal-to-noise eigenvalues allow also a compression of the data 
(see E50]). 

In addition to the CMB temperature, the determination of the statistical 
nature of polarization is crucial to test the Gaussianity of the CMB. Most 
of the tests proposed above for the temperature can be implemented for the 
CMB polarization. There are, however, other properties which are specific of 
the pseudo- vector nature of polarization. Up to now only a few tests have been 
developed based on geometric characterisitics of polarization: peculiar points, 
up-crossing and down-crossing points of the modulus of polarization (see for 
instance |171| ). In any case, it would be very desirable to device methods 
combining temperature and polarization maps in an optimal way. 

Some of the methods already discussed can as well be used to search for 
possible features present in the CMB maps arising from global anisotropy, in- 
homogeneity or topology. Particularly interesting are those based on wavelets, 
optimal filters and pattern recognition. If the universe is anisotropic/inhomogc- 
neous (as in the case of, e.g., Bianchi, Tolman or Swiss-Cheese models) then 
the CMB temperature fluctuation field will be anisotropic/inhomogeneous on 
the sphere. The space-time geometry of the universe is expected to affect 
mainly the low-order multipoles (i.e. large angular scales). However, focusing 
of geodesies can also produce a large variety of features at smaller scales (e.g., 
hot spots, rings of fire, spiral patterns, North-South asymmetries, etc., see 
|19l I1U9| . For recent applications to detect those features in WMAP data see 
Sec. 16.21) . If the topology of the universe is not trivial then anisotropic fea- 
tures will naturally appear in the CMB maps. Also symmetries are expected in 
some topologies (e.g. toroidal one). A method to probe non- trivial topology is 
based on the S-statistic which measures the symmetry of a map |228II66| . An- 
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Fig. 11. Spherical Mexican Hat Wavelet coeficients for COBE-DMR (from [IT]. 

other technique consists on the identified circles principle. Circular patterns of 
temperature fluctuations are generated in compact topologies. A given circle 
shoud be seen in different positions of the CMB sky since it can travel to the 
observer following different geodesies. Notice, however, that it will not look 
exactly the same because of the ISW effect caused by different intervining 
structures and the Doppler effect at recombination whose projection on dif- 
ferent geodesies will be also different. This technique has been applied to the 
COBE-DMR data Techniques based on wavelets, that allow to detect 

the patterns of spots imprinted by the different topologies, have been imple- 
mented in |197| . The results show the potential of scale-scale correlations for 
distinguishing among different topologies. Recently, other techniques involv- 
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ing phase correlations have been applied to the WMAP data to constrain the 
topology of the universe [7J . 

6 Data analysis beyond the power spectrum 
6.1 The Search for non-Gaussianity before WMAP 

The first systematic analyses that have been done on CMB data in search 
for non-Gaussianity, have been performed on the 4 year COBE-DMR data 
set. This was for a decade the only whole-sky data set publicly available. 
The COBE-DMR data have been found to be compatible with the Gaussian 
hypothesis by the majority of the statistical tests applied. The few analyses 
finding deviations from normality have been later proved to be either incom- 
plete or interpreted as undocumented systematic effects (see below). Monte 
Carlo simulations, taking into account all the instrumental and observational 
constraints of the data under analysis, arc usually performed to estimate dis- 
tributions of the testing-statistic as well as confidence levels. 

The degree of asymmetry of the data was measured with the 3-point corre- 
lation function (equilateral and pseudo-collapsed) by |12U| . Extrema distribu- 
tion and correlations of the temperature field have been tested by [1201 1183j . 
Several works have relied on morphological characteristics of the observed 
field to detect non-Gaussianity, the Minkowski junctionals, |2()7I I175| and in 
particular the genus (number of peaks above the threshold minus the number 
of holes. H2D1 CH31) to study the COBE-DMR data. The Partition Function, 
based on the combination of information at several scales and the contribution 
of different moments of the measure (defined in this case as the sum of abso- 
lute temperatures in a defined box), and the roughness of the last-scattering 
surface were proposed by [7Dj and |165| . respectively, as powerful statistics to 
detect deviations from Gaussianity. All the analyses performed in those works 
were done in real space. No evidence of departure from Gaussianity of the 
COBE-DMR anisotropies was found in all the cases discussed above. 

Tests were also carried out in spaces defined by eigenmodes and wavelets 
(see previous section). A Principal Component Analysis was done by [23 to 
extract the eigenmodes from the COBE-DMR data. \ 2 > Kolmogorov- Smirnoff, 
cumulants up to fourth order and significance of the top few outliers were cal- 
culated from the eigenmodes and tested against Gaussianity. None of these 
tests rejected the Gaussian hypothesis at the 95% confidence level. Skewness, 
kurtosis and scale-scale correlations of the COBE-DMR Wavelet coefficients 
have been studied in the north and south faces of the quad-cube COBE pix- 
clization, using different wavelet bases |1771 11691 15], Spherical wavelets have 
also been convolved with the COBE data (HEALPix pixelization) to test the 
distributions of wavelet coefficients at several scales against Monte Carlo sim- 
ulations assuming Gaussianity 1411 14ri| . A detection of non-Gaussianity 
was claimed by |177| finding the value of the scale-scale correlation (between 
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scales 11° and 22°) for the COBE data outside the 99% confidence level. 
However, |I69| pointed out the orientation dependence of orthogonal wavelet 
basis. A rotation of the data by 180° makes the COBE data set compatible 
with Gaussianity. 

Several works have studied the bispectrum of the COBE-DMR data 
in a simplified form that ignores correlations between different multipoles 
123 ESI El- A rejection of Gaussianity was claimed by jJ7| at the 98% confi- 
dence level and appeared at I = 16. In a subsequent paper J2j analysed the 
bispectrum as a function of frequency, channel and time interval. The non- 
Gaussian signal was found to come only from the 53 GHz map, implying a 
non-cosmological origin. Removal of data corresponding to a period of time 
dominated by a systematic effect (data collected during an eclipse) resulted 
in a COBE-DMR bispectrum compatible with Gaussianity. 
An extended bispectrum including components with correlations among differ- 
ent multipoles was studied by |147) . Gaussianity was rejected by the COBE- 
DMR data at a confidence level larger than 98%. An exhaustive study of the 
possible non-cosmological origin of this detection is included in that work giv- 
ing negative results. With the aid of the WMAP data, this detection was later 
proved to be due to pixelization artifacts of the Quad-Cube teselation of the 
sky and the eclipse data. 

All the previous works calculated the COBE-DMR bispectrum for a limited 
number of modes. More recently, |124| have measured all the modes available 
from those data, 466 in total. They conclude that the data are consistent with 
Gaussianity. A recent study using the trispectrum also shows consistency with 
Gaussianity (having previously removed the data affected by the systematic 
effect, [US]. 

While all the works mentioned above performed blind analyses, some works 
tried to detect deviations from Gaussianity in the form given by the non-linear 
parameter Jnl defined in Sec. 14.1 1 Only upper limits were obtained at least 
two orders of magnitude above the values expected in the most optimistic 
non-standard inflationary models (see e.g. |124l I43| for analyses involving the 
bispectrum and spherical wavelets, respectively). 

Uncertainty remained about whether the CMB anisotropics are Gaussian 
distributed or not. Scales observed by the COBE-DMR experiment are above 
the causal horizon at the last-scattering surface. It could well be that even 
if the anisotropics are non-Gaussian distributed at smaller scales, they ap- 
pear Gaussian above degree scales due to the central limit theorem. Therefore 
small angular scale observations are needed to answer this question. At present 
many of these new small scale observations have been tested against Gaus- 
sianity. The genus of the data of the balloon-borne QMAP, ground-based 
Saskatoon and QMASK (a combination of QMAP and Saskatoon) experi- 
ments was compared to Gaussian predictions by |178| . Only in the case of 
QMASK was an asymmetry in the genus curve found. However this devi- 
ation from the Gaussian genus curve was not statistically significant. |2()9j 
calculated several morphological quantities from the QMASK data and con- 
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eluded that they were consistent with the values expected for a Gaussian 
field. Several statistics have been computed to test the Gaussianity of the 
balloon-borne MAXIMA- 1 data ESB EEB El E| ■ HSU apply moments, cumu- 
lants, the Kolmogorov-Smirnof test, the x 2 test, and Minkowski functionals 
in eigen, real, Wiener-filtered and signal-whitened spaces; |203j estimate the 
bispectrum; @1 E] apply goodness-of-ht tests to the data. Consistency with 
Gaussianity was found in all cases. The balloon observations of BOOMERanG 
maps were also searched for deviations from Gaussianity. |188| perform a real- 
space analysis, computing skewness, kurtosis and Minkowski functionals of 
those data. No significant deviation from Gaussianity was found. 

More recently and contemporary to the WMAP analyses, the statisti- 
cal distribution of interferometric observations taken with the VSA experi- 
ment have been also examined. |205j apply some tests in real and harmonic 
spaces, the former to the maximum-entropy reconstruction of the observed 
regions. |218| focus on the bispectrum of the data. |2U1| apply smooth tests 
of goodness-of-ht developed by |7j . The results of the analyses seem to be in 
relatively good agreement with the Gaussian assumption (except for maybe 
one or two fields in the analysis of [201 ). 

6.2 WMAP and non-Gaussianity 




-0.44 ^^^^^^^^^^^^Hi i^^^^^^^^^^m <-l}.33 

Fig. 12. The first year WMAP combined, foreground cleaned Q-V-W map together 
with the kpO mask. 

The quality of the first-year all-sky data provided by the WMAP satellite 
has motivated many works dealing with the Gaussianity and isotropy of the 
CMB signal (also in relation to the geometry and topology of the universe, 
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see Sec. I4.2fl . Clean maps are constructed by the WMAP team and made 
publicly available to the community, together with noise templates and beam 
responses for all the channels. Also a mask is provided to ignore those pixels 
in the analysis where the foreground contamination is believed to be unavoid- 
able. All this information is essential to perform "realistic" simulations of the 
data needed to compare the values of the chosen statistics with the ones ob- 
tained from the true data. Fig. 1 121 shows the first-year WMAP map obtained 
after a noise-weighted combination of the foreground cleaned Q-V-W chan- 
nels (hereafter WCM). Although other foreground cleaned maps have been 
produced based on different techniques [2^1 12251 1751 1180| . for non-Gaussianity 
studies the WCM (together with the most conservative kpO mask) is preffered 
since the noise properties are well defined for this map. 

There have been a wide range of analyses of the temperature distribu- 
tion using different quantities. Minkowski functionals have been computed by 
[1251 1511 11791 11481 174| . The results obtained from those works are diverse: 
whereas compatibility with Gaussianity is found in |1251 l5*T| , others find sig- 
nificant deviations from it (1751 174] . The main reason for the discrepancy 
is that, although based on the same quantities, they construct different mea- 
sures from them as the asymmetry in the genus above/below a given threshold 
or in different hemispheres (these latter asymmetries have also been found in 
the power spectrum or the n-point correlation functions for different regions 
of the sky, see below) . 

N-point correlation functions have been calculated by (HHUZH] for the tempera- 
ture. jZH] find North-South asymmetries on large scales in ecliptic coordinates. 
Statistics based on the alignment of low-multipoles have been computed by 
HEIES1E221EIZIE1SIII301ES1 Unconfortably low probability for the observed 
quadrupole and octopole alignment has been obtained in those works (the so- 
called axis of evil) . The bispectrum have been calculated for different subsets 
of all its components up to the WMAP resolution [1251 11321 135] . Significant 
non-null results have been obtained in several of those analyses which agree 
well with the emission expected from undetected point sources [1251 ITU] . In 
addition those results can be transtated to limits on the non-linear parameter 
Jnl [1251135] . By using another subset of multipolc configurations, namely the 
closest inter-^ components, 132 have also found North-South asymmetries in 
agreement with previous results (7J3 E2J • 

Phase correlations have been used to test the hypothesis that the CMB 
temperature fluctuations on the celestial sphere constitute a homogeneous and 
isotropic Gaussian random field. This hypothesis implies that the phases of 
the spherical harmonic coeficients corresponding to the temperature fluctua- 
tions are independent and uniformly distributed in the interval [0, 2w]. Phase 
analyses of different derived maps of the WMAP data have been performed 
by [49111721 1501I173| . These works find departures from uniformity and/or in- 
dependency which in general reflect galactic contamination or noise gradients 
rather than primordial non-Gaussianity. 
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Several statistics based on the spherical wavelet coefficients of the CMB 
data have been used to test the Gaussian hypothesis [2381 117U1 1161)1 IBlfl FTTl 
1143] . Highly significant deviations have been found in the kurtosis of the 
WMAP data for both symmetric [238U17U) and directional wavelets |16(Jj . An 
analysis of extrema of the wavelet coehcients identifies a very cold spot as the 
cause of the large kurtosis found with the spherical Mexican Hat wavelet [S2| 
(see Fig. ll3[) . This cold spot has been also revealed by the higher criticism test 
applied to the wavelet coefficients 03]. Moreover, j^S! find a nearly spherical 
shape for the cold spot and also that it cannot be explained with foreground 
emissions. A test of global isotropy has been recently performed by |248| using 
steerable wavelets. Convolving the CMB map with the steerable wavelet given 
by the second derivative of the Gaussian at a fixed scale, a prefered orientation 
is signaled by the maximum of the coefficients corresponding to the different 
orientations and assigned to each pixel. Each pixel is thus looking at all the 
pixels intersected by the great circle formed by extrapolating the prefered 
orientation. Then the number of times each pixel is seen by any other pixel is 
counted. In the case of a homogeneous and isotropic Gaussian random field 
all the pixels are seen the same number of times on average. Performing this 
test on the WMAP first-year data, a number of anomalous directions were 
found at a very high significance level. Their mean direction is very close to 
the ecliptic poles, being w 10° from it (see Fig. ll4f> . Moreover, the anomalous 
directions are located along a great circle whose perpendicular is pointing 
towards a direction close to the dipole and to the axis of evil. This result 
synthesisc for the first time the previously reported anomalies commented 
above on the North-South asymmetry and the alignment of low-multipoles. 
Contrary to these results, analyses based on the bipolar spectra do not find 
violation of the global isotropy of the universe EH ■ 

Number and correlation of extrema are computed by jl84lllH5ll2^2| . [1341 
1135) find that the maxima and minima in the WMAP data are not hot and 
cold enough. 232 find evidence for non-Gaussianity on large scales which 
could be associated with residual foregrounds. 

Local curvature has been investigated by |93l l3"5| . They classify the points 
in the data as hills, lakes or saddles according to the sign of the eigenvalues of 
the Hessian matrix, and consider the proportions of each of them above a cer- 
tain threshold. A clear asymmetry is found which is maximum for hemispheres 
centered near the ecliptic poles, consistent with the results by [75], 

Several suggestions have been made to explain the anomalies found in the 
first-year WMAP data at the large angular scales, namely the power deficit 
and alignment of the low-order multipoles, the asymmetries between different 
hemispheres and the cold spot 4 . 

A suggestion to explain the low-multipolc anomaly invokes a non-trivial 
topology which due to the finite size of the universe tends to cut power at the 
lowest multipoles and in some cases to align them [1441 11991 ITT] , However, 



4 These anomalies remained almost unaltered in the 3-years WMAP data. 
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WCM, scale = 300' 




Fig. 13. The Cold Spot found in the WMAP combined map after convolving with 
a Spherical Mexican Hat Wavelet with scale 300' (see |238M62M63p 




Fig. 14. Excess of the number of times each pixel is seen by any other in a units, 
where a is the dispersion corresponding to Gaussian simulations (see |248| 'l 
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apart from this nice consequence, there is no evidence for a possible non- 
trivial topology. Another possibility proposed is the lensing of the dipole by 
local large-scale structures which can potentially produce a quadrupole and 
octopole similar to those found in the WMAP data 236 . It remains to be 
demonstrated that the appropriate physical characteristics of those structures 
(in terms of mass, size and position) are in fact present in the local universe. 

Possible explanations for the cold spot are the RS effect produced by the 
most prominent voids or massive structures in the universe |151l 12331 ll()7j , 
topological defects in the form or textures |234| or, a more radical one, based 
on an inhomogeneous model of the universe given by a finite sphere of dust 
and dark energy 

Anisotropic and inhomogeneous geometries leaving spiral or stripy pat- 
terns in the CMB sky have been proposed to correct the large-angle anoma- 
lies found (asymetries, the cold spot and the low-multipole power deficit and 
alignment) [1081 11311 1451 132) (see however |162| for further non-Gaussian de- 
tections in the Bianchi Vllh corrected WMAP first year data). Although a 
good fit can be found for most of the anomalies with the anisotropic geome- 
tries (Bianchi Vllh) however the values required for the matter content are 

too high unniES- 

Small changes in the magnitude of the WMAP dipole vector |50). together 
with possible systematics or foreground residuals associated with the ecliptic 
(as might be infered from the recent results by |75l 1921 12*4"%] ~1. might be consid- 
ered as a possibility to explain the large-scale anomalies. However much more 
work is still needed in order to find a convincing case for them. 

7 Conclusions 

Cosmology is passing through a period of abundance of new data which im- 
ply a significant improvement on the precission in the determination of the 
parameters characterizing the universe. The CMB field is probably the major 
contributor to that improvement and is expected to be so in the next decade. 
With only the first year of WMAP data , the power spectrum of the temper- 
ature anisotropies has been determined with error bars below cosmic variance 
up to multipoles £ ~ 600. Information on larger multipolcs is already provided 
by several high-resolution experiments allowing a consistent determination of 
the third acoustic oscillation and providing an incomplete coverage of the 
multipoles up to £ 2500. In the next few years, the ESA Planck mission is 
expected to provide accurate and high resolution all-sky images of the CMB 
temperature anisotropies in a wide frequency range which will allow an op- 
timal removal of foregrounds and a very good control of systematics. From 
those high quality images a complete coverage of multipoles at the cosmic 
variance limit will be obtained up to £ 2500. The Cg so obtained will be the 
ultimate temperature power spectrum of the primary CMB anisotropies since 
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secondary anisotropies are expected to start to be dominant beyond those 
multipoles. 

The large-scale galaxy surveys as well as the high-z SN la light curves 
have also played a crucial role in establishing the concordance model. The 
strong complementarity of these two cosmological data sets with the CMB 
anisotropies has implied a significant improvement on the parameter determi- 
nation compared to each data set considered individually. And, what is more 
important, a consistency check of the best fit model. In the next years, that 
combination of data sets is expected to be strengthened with weak lensing 
observations of background galaxies in large areas of the sky. 

In addition to the power spectrum, the recent studies on the statistical 
distribution of the temperature anisotropies provide independent tests of the 
standard inflationary model. Significant deviations from a homogeneous and 
isotropic Gaussian random field have already been found in the first-year 
WMAP data. However the origin of those deviations, whether cosmological, 
residual foregrounds or systematics, is still a matter of debate. More precised 
maps from next years of WMAP data and specially the high quality maps 
expected from the Planck mission will certainly help in the clarification of the 
possible causes of the anomalies found. 

Finally, high sensitive polarization maps to search for the B-mode will be 
the next major quest after Planck. The B-mode polarization is at present the 
best way to study the primordial background of gravitational waves expected 
in the standard model of inflation. It is therefore not surprising that both 
space agencies, NASA and ESA, have considered polarization missions in their 
future programs "Beyond Einstein" and "Cosmic Vision" , respectively. 
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